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"Tunneling" Amplitudes of a Massless Quantum Field. 
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r~| ■ Abstract 

^ I We propose a method for the approximate computation of the Green function of a 

. scalar massless field subjected to potential barriers of given size and shape in spacetime. 

^ \ The potential of the barriers has the form V{(j3) = 5,{4>^ — </'o)^i is very large and 00 very 

close to zero, the product ^^cf^ being finite and small. This is equivalent to the insertion 
of a suitable constraint in the functional integral for (j). The Green function contains a 
double Fourier transform of the characteristic function of the region where the potential 
has support. 

03.70.+k Theory of quantized fields. 
31. 15. Kb Path-integral methods. 

Several problems in Quantum Field Theory require the introduction of constraints or bound- 
ary conditions on the field. The best known case is the Casimir effect [1]. 
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Usually problems of this kind are reduced to partial differential equations for the propagator 
of the field (and thus to eigenvalues problems), with suitable boundary conditions [2]. We would 
like to show that in certain cases it is possible to obtain, to lowest order at least, an explicit 
solution, that can be evaluated numerically. 

Let us consider a massless scalar field 0, coupled to an external source J. The system is 
described by an Euclidean functional integral, of the form 



W[J] = J exp |- J (fx {d(t)Y - ^'[0, J] j (1) 

Now we want to impose on the field the condition 

0^(a;) = 00 xeVti^ i = l...n, (2) 

where the fij's are regions in spacetime. In order to fix the ideas, we can suppose that the 
regions Qi represent compact, distinct portions of space, and are thus infinitely extended in 
the temporal direction. In an electrostatic analogy, they could represent perfect conductors of 
arbitrary shape, fixed at certain positions in space. 

In order to implement condition (0), let us choose the source term 

5'[0, J]=^J d^xJnix) [02(x) -(f)l]\ (3) 

where Jn{x) = 1 for xefli, Jn{x) = elsewhere. In the limit ^ — oo, the exponential in 
acts in the functional integration (|l|) like a (5-function and the square of the field, (j)'^{x), is 
forced to assume in each configuration the value 0q within the regions Qi. 

Expanding the square in (H) and inserting S' into (|l]) we obtain for W[J] the expression 

W[J] = J d[<l)] exp |- 1 d'x [id<pf - 2^<f)lJn{x)<l)\x) + ^Mx)<p\x) + ^Mx)4] } . (4) 

The last term in the square bracket is constant with respect to (pix), and its exponential can be 
factorized out of the functional integral. In a first instance - for weak fields - we can disregard 
the 0^(a;) term. We are then led to consider a quadratic functional integral, and the "modified 
propagator" G{x,y) = {(f){x)(f){y))j, which by definition satisfies the equation 

dl + jMx)] G{x, y) = -(27r) V(x - y), (5) 

where 7 = 2^0q > 0. Let us focus on the case when 0q = inside the regions f2j and let us take 
the limit 0o and ^ — > 00 in such a way that 7 is finite and very small, so that the term 



^yJaix) in eq. (j^) constitutes only a small perturbation, compared to the kinetic term. Then 
we can set 

Gix,y) = G\x,y)+^G'ix,y), (6) 

where G^{x,y) is the propagator of the free scalar field, and we find immediately that G'{x,y) 
satisfies the equation 

dlG'ix,y) = -Mx)G'ix,y). (7) 

In general G'{x,y), unlike G^{x,y), will not depend only on (x — y), because the source 
breaks the translation invariance of the system. It will then be necessary, in order to go to 
momentum space, to consider the Fourier transform of G"(x, y) with respect to both arguments. 
We define G'{p, k) and Jn{p) follows: 

G'{x,y) = j d^P j d^ke'P''e'^yG\p,k) (8) 

and 

Mx) = I d'^pe'P^Mp), G\x,y) = J d^k "" . (9) 

The right hand side of (0) can be rewritten as 

Mx)G\x,y) = Jd'pJ d^ke^^JM'—^ = j '^^^j d'^pe'^ye'^'-^^^^ (10) 

and we obtain the following algebraic equation for the double Fourier transform of the first 
order correction to the propagator: 

pG{p,k) = — . (11) 

Transforming back, in conclusion we find 

G\x, y) = J d^p J d^k e^P^'e^'^y (12) 

Therefore, if we know the Fourier transform of the characteristic function Jq of the spacetime 
region where the constraint is imposed, we can in principle compute the leading order correction 
to the propagator of the field and thus to 

Now, it is known [3] that the vacuum-to-vacuum amplitude W[J] = (0^|0~)j of a field 
system in the presence of an external source J is related to the logarithm of the energy of the 
ground state of the system: 

Eo[J] = -T-HnW[J], (13) 



where the functional integral is supposed to be suitably normalized and the source vanishes 
outside the temporal interval [— T/2,+T/2], with T eventually approaching infinity. (We use 
units in which h = c = 1.) 

By imposing that the field vanishes in certain compact regions Qi of arbitrary shape and 
position, and taking at the end the derivative of W[J] with respect to their relative coordi- 
nates, one can find in principle the forces present, due to vacuum fiuctuations, in a system of 
generalized "uncharged conductors". 

A more interesting application of ([T2|) occurs in the case when the field 0(x) also interacts 
with static pointlike sources placed at xi,X2...X7v. Namely, let us add a further, linear 
coupling term 5*^ to the action of the system: 

N 

Sq= d^xQ{x)(f){x), with Q(a;) = ^g/3(x-Xj). (14) 

i=i 

The energy of the ground state of the system corresponds, up to a constant, to the static 
potential energy of the interaction of the sources through the field 0. It is obtained, as before, 
from the functional average of the interaction term, computed keeping the constraint into 
account: 

Eo[J,Q] = f/(xi,...,x^) = -T-Hn{exp{-SQ})j. (15) 

Expanding ([l5|), one finds that to leading order in the g^s, f/(xi, ...,XAr) is given by a sum of 
propagators integrated on time: 

[/(xi, ...,Xiv) = -T ^ Ijli J I dti{(t){tj,yij)(t){tu:i^i)) J (16) 

where tj,ti e [— T/2, +T/2]. Since the regions Vti are infinitely elongated in the temporal direc- 
tion, the function J^ij) + k) gets factorized as 

Mp + k) = {2ti)%p^ + A;o)L(p + k). (17) 

Clearly the potential is disturbed by the presence of the "barriers" ^^(x). We can write, to 
first order in 7, 

f/(xi,...,x^) = f/°(xi,...,x;v)+7f/'(xi,...,x^). (18) 
and taking into account eq.s (|T2l), ([T7|), we find 



f/'(xi, ...,XAr) = -T ^ qjqi I dtj I dtiG'{xj,Xi 



= -(2vr)^T-i f: q,q, j dt, j dU j d'p J d'k 6-^^-+-^^^' + k) ^ 

N I. I. I. I. pipo(tj-ti)+«Px,+«kxi;; / _l_ 

= -(27r)8r-i V mi / dt, / dti / d^p / d\^- — , , , .^ • 

Changing variables tot = tj—ti and s = tj+ti and integrating, we finally obtain the contribution 
from the perturbation to the static potential energy: 

t/'(xi, ...,x^) = -(27r)« J2 m JdpJ rfke^P-^+^'^''' (19) 

The six- dimensional integral (0), although formally simple, is in general impossible to 
evaluate analytically, because the Fourier transform of the constraint depends in a non 
trivial way on all angles, not just on |k| and |p|. In several cases, however, it is possible 
to rescale the integration variables dividing them by the sizes of the regions Qi and by the 
distances between the points Xj and x/. In this way, any power-like dependence of U' on the 
various lengths involved is factorized in front of the integral; the weak residual dependence 
(usually logarithmic) can be estimated by computing the rescaled integral numerically. 

In order to illustrate this point, let us consider now a special choice for j^. Clearly, if Jq 
was exactly equal to the characteristic function of the set Q (i.e., a step function in coordinates 
space like above: jn{'^) = 1 for x el], jn{'^) = otherwise), then the Fourier transform jn(p) 
would have a strongly oscillating behaviour. To make the integral (|T9|) more accessible to a 
numerical computation, it is useful to regularize the step of jn{'^) and "smooth" it. As a further 
simplification, let us consider a gaussian function: 

this represents a single region Q, centered at the origin, with approximate sizes a, b and c in 
the three directions. The transform of (EDI) is 



jn{p) = 7r3/2(a6c)exp 



1 / 2 2 1,2 2 2 S 

j[-a p^-bpy-c p^ 
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The generalization to the case of several regions Qi, centered at the positions Xj, is obvious: 

■1 



(21) 



jf,(p) = ^7r3/2(a,6,Q)e^P^'exp 
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-a Px ~ b py - c p 



(22) 



After inserting (^) into eq. (p!9l), the integration variables can be rescaled in such a way 
that the integral contains only the ratios between a, b and c. The dependence of the integral on 



these ratios is quite weak; the main dependence is given by the factor (abc) and by the factors 
coming from the re-definition of the integration variables. More detailed numerical results will 
be presented elsewhere. 

We observe that the representation of the constraint on the field through a "smoothed" 
jcii'x.) modifies in part the original physical interpretation of our formalism. Instead of perfect 
conductors, we are representing here more realistically some smooth potential barriers in the 
functional integral for the field 4>. The probability that the field survives, by tunneling, up over 
the barriers, depends on the parameter 7. 

We have seen that the constraint produces only a small perturbation in the interaction 
energy of pointlike sources. Thus our representation of the constraint by a potential barrier is 
suitable to describe non-perfect "conductors", or conductors whose size or thickness is much 
smaller than the wavelength characteristic of the interaction between the pointlike sources. 
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